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Abstract 

The Born approximation, one photon exchange, used for DIS is subject to virtual 
radiative corrections which are related to the long-range Coulomb forces. They may be 
Oh. sizeable for heavy nuclei since Za is not a small parameter. So far these corrections are 

r- 1 \ known only for two processes, elastic scattering and bremsstrahlung on the Coulomb 

field of a point-like target. While the former amplitude acquires only a phase, in 



X 
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the latter case the cross section is modified also. Although the problem of Coulomb 
corrections for DIS on nuclei is extremely difficult, it should be challenged rather than 
'swept under the carpet'. The importance of these radiative corrections is questioned 
in present paper. We show that in the simplest case of a constant hadronic current 
the Coulomb corrections provide a phase to the Born amplitude, therefore the cross 
section remains the same. Inclusion of more realistic hadronic dynamics changes this 
conclusion. The example of coherent production of vector mesons off nuclei reveals 
large effects. So far a little progress has been made deriving lepton wave functions 
in the Coulomb field of an extended target. Employing available results based on the 
first-order approximation in Za, we conclude that the Coulomb corrections are still 
important for heavy nuclei. We also consider an alternative approach for extended 
nuclear targets, the eikonal approximation, which we demonstrate to reproduce the 
known exact results for Coulomb corrections. Calculating electroproduction of vector 
mesons we again arrive at a large deviation from the Born approximation. We conclude 
that one should accept with caution the experimental results for nuclear effects in DIS 
based on analyses done in the Born approximation. 
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1 Introduction 



Smallness of the fine structure constant a usually justifies lowest order perturbative QED 
calculations. In some cases, however, the expansion parameter is not small, for instance, for 
interactions with heavy nuclei, where Za ~ 1. Therefore, the validity of Born approximation 
for deep-inelastic scattering (DIS) should be questioned since the incoming and outgoing 
leptons propagate in the Coulomb field of the target, as is illustrated in Fig. [I], which can 
cause a deviation of the DIS cross section from the Born form. This may lead to important 




Figure 1 : Deep-inelastic electron-nucleus scattering. The thin wavy lines 
illustrate the long-range Coulomb forces which violate the DIS kinematics 
(q Pi ~ Vi) aD d modify the lepton wave functions. 

modifications of experimental data for nuclear effects in DIS which rely on analyses based 
on the Born approximation. 

The simplest example is the elastic lepton scattering in the Coulomb field of a nucleus. 
Neglecting the nuclear structure, i.e. treating the target as a point-like charge, one can 
solve this problem exactly. Gordon and Mott |2| have done it in the framework of 
nonrelativistic quantum mechanics. It turns out that the resulting amplitude is different 
from the Born approximation only by a phase which has opposite signs for positive and 
negative leptons. Thus, the elastic cross section has no Coulomb corrections. On the other 
hand, in the case of hadronic elastic scattering the Coulomb effects do modify the cross 
section because of the relative Coulomb phase || f|. [| . 

In the relativistic case the problem of elastic Coulomb scattering was solved for spinless 
particles by Schrodinger || and for fermions by Mott and Darwin [^]. Although the 
scattering phases were calculated, summation of the partial amplitudes, i.e. construction of 
the wave function of a relativistic charged particle in the Coulomb field, is still a challenge. 

Nevertheless, for practical applications one can employ the approximation of high orbital 
momenta, L 2 ^> (Za) 2 , which helps to solve the problem. Corresponding wave functions in 
the Coulomb field were found by Furry || and Maue and Sommerfeld ||. 

Bethe and Maximon ]10|] made use of the same approximation to solve a more compli- 
cated problem of bremsstrahlung and production of electron-positron pairs by an electron 
scattering in the Coulomb field. Their finding is interesting, the long-range Coulomb forces 
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modify not only the phase, but lead to a suppression of the cross section compared to the 
Born approximation. The effect is significant, about 10 — 15 % for heavy nuclei. The fact 
that the Coulomb forces suppress radiation can be intuitively understood as a manifestation 
of the Landau-Pomeranchuk effect. Indeed, the electron trajectory is bending smoothly over 
long distances, and the photons radiated from different parts of the electron path interfere 
destructively. 

One may expect similar modifications in other reactions, especially those which involve 
strong interactions, like DIS, they should be revisited aiming to clarify the role of long- 
range Coulomb forces. All analyses of DIS data are based on the Born approximation for 
the electromagnetic interaction with the target, i.e. the wave functions of the upcoming 
and outgoing leptons are treated as plane waves. At the same time, precise measurements 



TTJ disclose rather fine, few percent nuclear effects, like small antishadowing at x ~ 0.1, or 
variation of nuclear shadowing with Q 2 . It is still questionable how these small effects are 
related to the missed distortions generated by the Coulomb field. 



On the other hand, the HERMES experiment [E| has discovered recently unusually 
strong nuclear effects which cannot be understood within conventional approaches . 
Since these effects arise only after the radiative corrections are introduced into the analysis, 
it also motivated us to have a closer look at the problem of reliability of radiative corrections. 

This paper is organized as follows. In Sect. [| we start with a specific simplified example 
of DIS with a hadronic current independent of Q 2 . In this case the amplitude is proven to 
gain only a phase leaving the cross section unmodified. This is a new nontrivial result which 
we utilize for further applications in this paper. 

As a more realistic, but still simple example we choose the coherent vector meson pro- 
duction off nuclei which is treated within the vector dominance model (VDM). In Sect. fJTI 



we consider the case of a point-like target and arrive at a substantially modified cross section 
which deviates from the Born approximation. The Coulomb corrections increase with Q 2 
as is shown in Fig. |2|. 

The next step toward a realistic hadronic current is done in Sect. 3^2. It corresponds 
to the vector meson production off an extended nucleus, but with the same leptonic wave 
functions calculated for a point-like Coulomb center. Again large corrections to the cross 
section are found whose value and sign varies with lepton energy and Q 2 , as is demonstrated 
in Fig. [3|. 

Eventually, in Sect. [0]we replace the exact wave functions of the lepton in the Coulomb 
field of a point-like charge by approximate ones calculated for the case of an extended nuclear 
target. Unfortunately the only solution available in the literature corresponds to the first 
order expansion in Za. Assuming that it is not a phase we arrive at a significant correction 
up to 13% for heavy nuclei. 

In Sect. f|try a different technique, the eikonal approximation, which is designed in a way 
which allows to perform calculations for nuclei of large size. First of all, we check whether the 
eikonal approximation is able to reproduce the Coulomb corrections previously calculated 
exactly. We demonstrate in [Appendix B| for the case of e + e~ production off nuclei that the 
exact results for the Coulomb corrections are fully reproduced by the eikonal approximation. 
We calculate the Coulomb corrections for the process of vector meson production in the 
eikonal approximation in the limit of very high energies. Again the corrections are found to 
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be large and to saturate as function of Q 2 at Q 2 > 0.1 GeV 2 , as is shown in Fig. |[ 

Although our results are too rough to be incorporated into analyses of data for DIS, 
we conclude in Sect. [| that one should accept with a precaution the available experimental 
results for nuclear effect in DIS. We have found only one specific case of an oversimplified 
hadronic current when the Coulomb corrections have the form of a phase. In more realistic 
situations the hadronic current consists of a few (or many) terms which acquire different 
phases and lead to a modified cross section. Thus, all available data for DIS on heavy nuclei 
based on analyses utilizing the Born approximation might change if the Coulomb corrections 
are applied. 



2 A little theorem 

This section is aimed to demonstrate that for a specific simple form of the conserved hadronic 
current the effects of long-range Coulomb forces are reduced to a phase factor, like it happens 
for elastic scattering. This result then will be implemented into more complicated situations. 

Let us start with a process which leaves the nucleus intact, for example coherent elec- 
troproduction of vector mesons. The amplitude of this reaction in Born approximation 
reads, 

MieA^e'VA)^ [ d 4 xd*y (p 2 \j^x)\ Pl ) D, u (x - y) (V|J M (y)|0) , (1) 



where 

(p 2 \j„(x)\pi) = % 2 (x)'Y»% 1 (x) (2) 

is the operator of lepton current; J^iy) is the hadronic current operator; D^ v {x — y) is the 
Green function of the photon. Using the Feynman gauge we have, 

„ e iQ{x-y) 

D^tx -y) = -g^ u J d A Q ^ 2 - — , (3) 

and 

M(eA^e'VA)^-J ^9_j^Q ]Pu p 2 ) J„(Q ; P v ) . (4) 
Here the Fourier transform of the currents reads, 

jn(Q\Pi,P2) = J d 4 x(p 2 \j^x)\ Pl )e iQx ; (5) 
UQ;Pv) = Jd 4 y(Pv\Uy)\0)e- iQ »; (6) 

Pi,2 are the initial and final lepton momenta; Py is the momentum of the produced V . 
If the initial and final wave functions of the lepton are undistorted plane waves, 

V Pl , P2 (x) = e tpi < 2X u( Pl , 2 ) , (7) 

then 

3n(Q',Pi,P2) = (2n) i 5(p 1 -p 2 -Q)u(p 2 )-f tl u(p 1 ) , (8) 
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and 



M(eA -> eti) oc 



■ u(j> 2 h^u(pi) Jpipz - Pi; Pv) 



(9) 



(pi - P2) 2 

However, in the presence of a Coulomb field one must rely on the solution of the Dirac 
equation, rather than on the plane waves, 

. d 



—id-V + /3m + e> 



(10) 



Here a and (3 are the standard Dirac matrixes; 0(r) is the Coulomb potential which is 
independent of energy in the rest frame of the nucleus; m is the lepton mass. The initial 
state leptonic wave function which is a solution of this equation, has the form of sum of 
a plane and a spherical outgoing waves, \I/ + (pi,f) at r — > 00, while the final state wave 
function should contain a plane and an incoming waves, ^~{p 2 ,r). 

In general case equation ([TUP can be solved only numerically, but if the effects of nuclear 
size can be neglected, i.e. 

Ze 



an approximate analytical solution exists, as it was found by W.H. Furry 



ty~(p 2 ,r) 



c 



C* 



,W2-r 



1 - 



id ■ V 
id ■ V 



F [iZa, 1; i{p x r - pi ■ r)] u(pi) , 
F [—iZa, 1; -i(p 2 r + p 2 • r)] u(p 2 ) 



(11) 

(12) 
(13) 



where 



C 

e l,2 



e nZa/2 T(l-iZa) 



p? 2 + m 2 , 



(14) 



u{pi, 2 ) are the 4-spinors of the leptons, Y{x) is the gamma-function, and F(a, b; c) is the 
confluent hypergeometric function. The condition that the leptons are ultrarelativistic, 
ei,2 3> in, implies that the essential orbital momenta are large, / ^> 1. It is well satisfied in 
all cases we are interested in. 

For the sake of simplicity we will work with "spinless leptons" since we are interested 
only in an estimate of the effects and do not expect a principal modification related to the 
lepton spin. Then the lepton wave functions satisfy the Klein-Gordon equation, 



A*(r) = {[e-e^)(r)f-m 2 } f(r) 

= p 2 - 2 e e (f>(r) + e 2 2 (r)| #(r) . 



(15) 



Apparently, at high energies the term e 2 2 (r) can be neglected, and for the Coulomb po- 
tential Eq.flTT|) one can get the exact solution, 



* + (pi,r) = CF[zZa,l;z(pir-pVr)]e^- f 
ty~(p 2 ,r) = C* F [-iZa,l; -i(p 2 r + p 2 ■ r)] e 



(16) 
(17) 
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Correspondingly, for the lepton current, j^x) = j^(r,t), 



J(e 1 -e 2 )t 



where 



jo{r) = e^-*{p 2 ,r) (e x + e 2 + ^ + (Pi 
j(r) = ie 



This current is conserved, 



V^" (p 2 ,r) ^ + (pi,f) - (p 2 ,f)V^ + (p 1 ,f] 
dj (r,t) 



dt 



or, in momentum representation, 



vjo{Q) =q-j(q) 



(18) 

(19) 
(20) 

(21) 
(22) 



where u, q are the energy and momentum transferred to the target, v 2 — q 2 = —Q 2 , and 
Fourier transform of the current reads, 



(23) 



j IM (Q) = Jd 4 xe^j,{x) . 
In terms of time-independent current it reads, 

j (Q) = 2ir5{e 1 -e 2 -u)Jd 3 re i * r j (r) 



j(Q) = 27r5( ei -e 2 -u)J d 3 re l ^j(r) . 
Using these expressions and conservation of the hadronic current, 

vJ Q {Q) = q-J(Q) , 
we arrive at the following form of the amplitude Eq. @, 



M = 2vr 



d 3 q 



-Q 2 - io 



v 



J{Q) 



(24) 
(25) 

(26) 
(27) 



The integral in Eq. (|27|) contains hadronic current J(Q) with unknown dependence on 
Q [Q dependence of the leptonic current is fixed by Eqs. (^4|) - ( p5j) 1. As the first simple 
trial we assume that it is Q independent, J{Q) = c. Then the amplitude Eq. ( |2~7l) takes the 
form, 

M = 2i (2vr) 2 — , (28) 



where 



1 



d 3 r 



(29) 
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Let us consider kinematics when the virtual photon takes a finite fraction y = v/e\ of 
the initial lepton energy, and the scattering angle is small, 9 <C 1. In this case the relative 
contribution of the the term Za/r in the current Eq. ( |19"|) estimated with the plane wave 
approximation is small, ln(l/6>) 9 2 /y <C 1. Thus, it is suppressed by a factor m^/V. 
Neglecting this term we simplify the current Eq. fll~9|) to 

3o(r) = (ei + e 2 ) f) #+(pi, r) . (30) 

Applying Eqs. (|T§ - (0) to Eqs. fl3§) - (|3Dp we get for the vector dEq. 



d = 
where 



d 3 r 



ivr+iq-r 



[em - e 2 pi) Fi(f) F 2 (r) + i Cl ^(f) VF 2 (f) + i e 2 Viq(r) F 2 (f) 



Fi{r,Pi) = CF[iZa, l;i(pir -pi • r)] , 
F 2 (f,p 2 ) = C* F[—iZa,l;—i(p 2 r + p 2 ■ r) 



Then we make use of the following relations, 



Pi 



V?F 2 {r,p 2 ) 



Pi 



Vp 2 F 2 (f,p 2 ) 



<iA e 



-Ar 



(31) 

(32) 

(33) 
(34) 

(35) 



and (see in fi"4|| ) 



d 3 



/(XT 
Fi(r,pi) F 2 (r,p 2 ) exp(iq ■ f - Ar + ivr) = I(q,p 1: p 2 , A) 
r 

F(iZa, 1 — zZa; 1; a;) , (36) 



where -F(a, 6; c; <i) is the conventional hypergeometric function, 

uv 



X 

w 
u 

V 

z 
N 



1 - 

wz 

f + (A - ivf ■ 

ipi + q) 2 -(P2 + v + ; 

-(Pi - q) 2 + (pi + is + i\) 2 ; 

(Pi + P2 + v + iA) 2 - (px -P2~q) 2 
2 irZa 



T(l-iZa) 



sinh(7rZa) 



(37) 
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With Eqs. 



we arrive at a new expression for d in Eq. ( j3l| ' 



d = (eip 2 - e 2 px) I(q,pi,p 2 , A = 0) 



q=pi ~P2 



- / d\g{p 1 ,p 2 ,X) 



(31 



(39) 



where 

g(pi,p 2 ,X) = i(p 2 eiVf 2 -p 1 e 2 'Vf 1 )I(q,pi,p 2 ,X) _ _ _ 

g=pi-P2 

Further calculations are moved to Appendix A where they are performed for a pho- 
ton of mass m 7 . In the final Eq. ( A. 26 ) the second term in the brackets vanishes since 
F(l + iZa, 1 — iZa\ 2;x ) diverges logarithmically, but 1 — x oc m 2 — > 0. The first term 
at m 7 — > gets, 



F(iZa, —iZa; 1; x 



m-y— »0 



|r(l-iZa)r AT 



(40) 



Then Eq. ( |A.26|) leads in the limit m 7 — > to 



- 4:71 



Q 2 



(Pi + £>2 + f ) 



(41) 



Thus, we have proven that the amplitude is different from the Born one only by the phase 
factor. 



3 Coherent electroproduction of vector mesons 
3.1 VDM, a point-like nucleus 

Apparently, the hadronic current depends on Q 2 . For instance, the VDM suggests 



ni 



J (Q) = AO) , ( 4 2) 

v 



Q 2 + m 



where my is the vector meson mass. Including the propagator of the virtual photon 1/Q 2 
it can be represented as, 

< = J_ l — (43) 

Q 2 (Q 2 + m 2 v ) Q 2 Q 2 + m 2 v { ) 

Correspondingly, the new the vector d equals to the difference between Eq. ( |3~8D and the 



same expression, but with replacement u => \Jv 2 — m 2 ?. As a result, on top of the phase 
factor the amplitude acquires another complex form factor, 

S ^ = 1 - - uZ \ xl ~ iZa W ( K ) > ( 44 ) 
smh(7rZa) 



S 



where 



W(k) = F(iZa,-iZa;l;K) - iZa(l - k) F(l + iZa,l - iZa;2; k) ; (45) 
Q 2 

Q 2 + m 2 v - 

We calculated the ratio of the cross sections of reaction I A —> I' V A calculated with the 
distorted Coulomb wave functions and in the Born approximation. The results are plotted 
it in Fig. as function of Za for different values of Q 2 /my. Deviation from unity increases 
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Figure 2: Deviation of the cross section from the Born approximation as 
function of Za for different fixed values of Q 2 /my. 



with Q 2 , and is of course larger for heavier nuclei. Note that in the limit Za — > the form 
factor Eq. ([0]) recovers the conventional VDM form Eq. (^) and the deviation from the 
Born approximation vanishes, as Fig. confirms. 

If the hadronic current contains a higher power of Q 2 , 

J <«> = J ^SFw)" " (46) 

it can be treated in a similar way provided that n is integer. In this case the previously 
obtained expressions should be differentiated in the parameter my. 

3.2 More realistic hadronic current 

We assumed above for the sake of possibility of analytical calculations that the hadron 
current J(Q) is independent of other variables, but Q. This means, in particular, that the 
amplitude of virtual coherent photoproduction 7* a — > V A is isotropic, i.e. independent of 
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the momentum transfer A = Q — Py. This is quite unrealistic suggestion, and a better form 
of the hadronic current would be 



J(Q,A) 



ey my 



1 



m 2 v + Q 2 1 + 5AV2 



(47) 



where ey is the polarization vector of the vector meson, and B is the slope of the transverse 
momentum distribution for the reaction — > V A. 

In this case all multiple integrations in the amplitude of electroproduction of vector 
meson can be done analytically down to the one dimensional integral, 



M(IA -> I'VA) = J dxe v ■ h(p u p 2 ,p v ,x) 
o 



where h — hi — h 2 , and 
hi, 2 (pi,P2,Pv,x) 



62 Pi V Pl - eip 2 V P2 + (eip 2 - e 2 pi 







dcoi,2 



q=pi-p2-xP v 



Here 



X = 


uv 
wz 




w = 


->2 2 

q -u ; 




u = 


(P2 + q) 2 ~ (P2 + uj) 2 ; 




V = 


(pi + uf + (p!-q) 2 ■ 




z = 


(P1+P2 + U) 2 - (pi -f 2 




UJi = 


(1 — X) V 2 — X(l — X) Py - 


-xB/2 ; 


UJ 2 = 


(1 — x){y 2 — m v ) — x(l - 


-x)p v -xB/2 



and 



J(pi,P2,q,Ui )2 ) 



d 3 r ■ 



(4J 



(49) 



(50) 



iujr+iq-r 



F(iZa, 1, ipiT — ipi ■ r) F(iZa, 1, ip2r + if>2 ■ r) 



4irN (w\ %Za 

— F(iZa, 1 — iZa; 1; x) . 

w V z J 



(51) 



We performed numerical calculations for forward production (A = 0) of transversely 
polarized vector mesons. The ratio of the calculated and Born cross sections cr/aBom = 
\M\ 2 /\MBom\ 2 is depicted in Fig. |3| for few values of ei at fixed y = (ei — e 2 )/ei = 0.6 and 
Q 2 = m 2 p . 

The magnitude of the effect is probably overestimated since the nucleus size is taken 
into account only in the production amplitude ^*A — > VA, while the wave functions of the 
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Figure 3: The same as in Fig. ^ but for different incident energies e x = 
200, 100, 30, 10, and 5 GeV at fixed y = 0.6 and Q 2 = m 2 p . 



incident and scattered leptons still correspond to the potential of a point-like center (see 
next Sect. |3.3| ). Nevertheless, these calculations reveal a new effect, the Coulomb corrections 
can be either positive, like in the previous case, or negative (compare with the results of 
calculations in 

We have also studied how these Coulomb corrections vary as function of y and Q 2 . 
While the magnitude of deviation from the Born cross section is nearly independent of y, it 
substantially increases with Q 2 , keeping about the same shape of the dependence on Za. 



3.3 Effects of nuclear size on the lepton wave functions 



One should expect the Coulomb effects to be diminished when the size of the nucleus is taken 
into account. Indeed, the V- meson interaction radius is short and the beam lepton has to 
have about the same impact parameter as the virtual photon or V-meson (the difference 
is ~ 1/Q). Propagating through the nucleus the lepton experiences a weaker electric field 
compared to the case of a point-like Coulomb center. Correspondingly, the Coulomb wave 
functions of the leptons must be corrected for the finite size of the nucleus. 

The general case of a massless lepton in an isotropic potential V(r) has been solved and 
the radial wave functions, correct in all orders in Za, have been found in [16] within the 



quasiclassical WKB approximation. They have been used later in [J17, [18|, to sum up the 
partial waves and build a formal expression for the lepton wave functions, 



^ ± (p,f) = e ±i&1 ' 2 n(r) e 



3 ±ib(J 2 -3/4) Jp-rr)(r) 



u(p) 



(52) 



where J = L + a /2 is the operator of the total angular momentum; 



2p 2 



dr 



1 dV(r) 
r dr 
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Ra 

S 1/2 = -Za ln(pR A ) - J drV(r) + b 



T](r) = — f dr' p — V(r) ; (53) 
pr J 1 J 
r o 

and u(p) is the 4-spinor of the lepton. Eq. fl52|) is a generalization of the Furry's wave 
functions Eq. (ji3|). 

Apparently, such an expression with the operator J in the exponent is not easy to use for 
practical applications. This is why only the first two terms of the order of (Za)° and (Za) 1 
in the expansion of the exponential in (|52"D have been considered in |17|, |l^, |19|. We skip 
here those lengthy expressions, but apply the procedure developed in [17], [TS| to our case. 



We assume a homogeneous charge distribution inside a sphere of nuclear radius Ra and use 
the approximation si^Ra 3> 1 and 6*12 <C 1, where 6 , i2 is the lepton scattering angle. Then 
the ratio of the amplitude modified by the Coulomb field to the Born one takes the form, 

M(IA -> VVA) . WZa _ . 3Za 

M Born (lA - I'VA) + l ~RT ~ ' ~5~ ' { ' 

Here 5 is the slope parameter of the differential cross section introduced in (|I7|). It is related 
to the mean charge nuclear radius squared, 

B= 1 -(r 2 ch } A = 1 -R A - (55) 
Thus, deviation from the Born cross section, 

a(lA^l'VA) 9 , 2 



1 + — (Za) 2 , (56) 



a Born (lA -> Z'VA) 25 

This is a sizeable correction for heavy nuclei, for example, the modified cross section on lead 
is 13% higher than the Born one. 

Surprisingly, the correction Eq. ([56]) does not expose any dependence on reaction kine- 
matics, on the contrary to the results of previous sections. It is probably a consequence 
of the higher order terms in Za missed in this calculations. How to elaborate with those 
terms is still a challenge which we leave for further studies. The purpose of the present 
estimates is to see whether the finiteness of nuclear size can substantially diminish the effect 
of long-range Coulomb forces. Apparently not, the deviation is still sizeable. One should 
calculate at least the next term 0(Z 2 a 2 ) to make it sure that the correction Eq. ( |5^ ) is not 
just a phase. We believe that such a possibility is very improbable, it happened so far only 
in the case of elastic scattering, and in the special artificial case when the hadronic current is 
a constant (Sect. |^). In both cases the Furry wave functions for a point-like Coulomb center 
are used. The example of Sect. |3TT] demonstrates that the even a simplest Q 2 dependence 
of the hadronic current leads to an amplitude which consists of few terms having different 
phases, and the cross section changes. 
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4 Eikonal approximation 



A different approach to the problem of distortion of the lepton wave functions in the Coulomb 
field of an extended nucleus is the eikonal approximation of Bjorken, Kogut and Soper (BKS) 
ipOj . It is not clear how precise this approximation is, the best way to figure it out is to 
compare its results with the known exact solution in the case when it is available. The 



exact cross section has been calculated by Davies, Bethe and Maximon (DBM) |21j for 
photoproduction of e + e~ pairs and Bethe and Maximon for bremsstrahlung ||10|| . Using the 
BKS approach we calculate in |Appendix B| the cross section for this reaction and arrive 
at exactly the same Coulomb correction to the Born cross section as it has been found in 
pT0|1 . This success can serve as an argument that the eikonal approximation is rather 
accurate. 

In the eikonal approximation the wave functions of spinless leptons read, 



exp 



ipi ■ r — i Xi\ r ) 



* (P2, r) = exp % p 2 ■ f + i X2(r) 



(57) 



Here the phase shifts are 



Xi{r) 



X2[r) 



dz' V(b, z') 



-oo 



dz'V(b,z') . 



(51 



At high energies, e 1)2 ^> y/Q^, the vectors pi, f>2 and pi — f>2 are nearly parallel. We chose 
the axis z along pi, correspondingly the vector f = (b,z) has projection b to the normal 
plane. 

For the long-range Coulomb field, V^(r) r ^oo = ±Za/r , the integrals in (|58|) are strictly 
speaking divergent. To fix the problem we introduce an infra-red cut-off, 



V(r) 



nr \ 



-fir 



where the last factor corresponds to the pole form of the nuclear formfactor with 

-2 _ ( r ch)A 
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(59) 



(60) 



One may interpret A in fl5"9"D as an effective photon mass, or (better justified) as the 
inverse screening radius of the nuclear Coulomb field by the atomic electrons. As soon 
as 1/A ^> Ra, any variation of the value of A may lead only to r-independent additive 
corrections to the phase shifts Xipi?) what does not affect the value of the cross section we 
are interested in. 

Utilizing the same realistic form of the hadronic current for electroproduction of vector 
mesons as in Sect. |3.2| we get the following expression for the amplitude, 



M(IA -»■ I'VA) = dxe v - f(puP2,Pv,x) , 



(61) 
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where 



/ — fi ~ h > 



and 



1 







2u±o dui.2 



eiVx 2 (r) + e 2 Vxi(^)]} 



x exp 



iff - ixi( r ) - i>X2\f) +iuJi 2 r 



(62) 



where k = f>\ — p 2 — xpy. 

Note that the ^-dependence of this expression comes via uj\^ defined in (|50|). The 
derivatives Vxi, 2 (?") in fl62| ) are the momenta transferred by the Coulomb field to the initial 
and final leptons. In the case of an extended nucleus they are of the order of Zch/Ra- 

Using relation 

(63) 



eij9 2 - e 2 pi 



eie 2 Q 2 



we see that for Q 2 3> R A 2 the second term in curly brackets in fl62"|) is much smaller than 
the first one and can be neglected. Then Eq. flolf) takes the form, 



(64) 



Xi(r")+X 2 (r)= J dzV(b,z) = x(b) = ±2Za{K (Xb)-K [^ + X)b]} 



where plus and minus correspond to different signs of the potential Eq. (|59|) . Since the sum 
Xi(r) + X2(^*) is independent of z one can explicitly perform integration over z in Eq. 
using the following relation, 



J — exp(i k l z + ito r) = 2 K (b \j k 2 l — lu 2 ^) , 



(65) 



where K (x) is the modified Bessel function, u = co>i )2 , and k^t are the longitudinal and 
transverse components of k. 

Thus, the problem of calculation of the amplitude (|61|) is reduced to integrals of the 
form, 

1 OO 

J dxj dbbK (b y/nl - uj 2 ) e ixih) J (M) , (66) 



which we calculate numerically. Our results for ratio of the calculated and Born cross 
sections are depicted in Fig. [| as function of Za. It turns out that the cross section is 
smaller than the Born one, and the difference is Q 2 independent at Q 2 > 1 GeV 2 . However, 
the range of Q 2 correlates with the photon energy v to satisfy the condition of coherence 
q L = (Q 2 + m 2 v )/2u<l/R A . 
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' 0.0 0.1 0.2 0.3 0.4 0.5 0.6 

aZ 

Figure 4: The same as in Fig. but calculated in the eikonal approxi- 
mation. The curve is calculated at Q 2 > 1 GeV 2 and is independent of 
Q 2 within the range of coherence qi < 1/Ra- 



5 Conclusions and outlook 

The Born approximation employed in all analyses of experimental data has no justification 
for heavy nuclei with Za ~ 1. Therefore, the accuracy of experimental results for nuclear 
effects in DIS may be essentially affected by this theoretical uncertainty. 

Radiative corrections are known to provide a substantial contribution to the DIS cross 
section. Usually they can be calculated rather accurately [E3], except one related to the 
modification of the lepton wave function by long-range Coulomb forces. This correction 
violating the Born approximation is usually ignored, and for a good reason: it is a very 
difficult task to calculate it. Nevertheless this problem should be challenged rather than 
'swept under the carpet'. 

A certain progress is made in the present paper towards the calculation of the DIS cross 
section using the lepton wave functions modified by the long-range Coulomb forces. We 
start with the simplest assumption that the hadronic current is independent of Q 2 and 
prove that the DIS amplitude acquires only a phase. In this case the Coulomb corrections 
do not modify the Born cross section. 

This conclusions, however, changes after a Q 2 dependence is introduced into the hadronic 
current. We demonstrate that on the example of coherent leptoproduction of vector mesons 
off a point-like nucleus which we treat within VDM. Of course VDM should not be used 
at high Q 2 where color transparency is important. However, we just want to estimate the 
scale of the correction, while precise calculations ready to use in an analysis of data is still 
a challenge. 

Further steps beyond the approximation of a point-like target lead to a conclusion that 
in the case of an extended nucleus the Coulomb corrections are still important and can 
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affect the existing experimental results for nuclear effects in DIS which are based on the 
Born approximation. To do more precise calculations one needs to know the lepton wave 
functions in the Coulomb field of an extended nucleus which are currently available only in 
the first order in Za. 

We also use the eikonal approximation which is better designed for the case of extended 
nuclei. First we demonstrate ( Appendix B ) that this approximation well reproduces the 



known exact results when they are available. Then we apply this method to the process of 
coherent electroproduction of vector mesons and again arrive at a sizeable correction to the 
Born cross section. 

We conclude that the long-range Coulomb forces may significantly modify the DIS cross 
section on heavy nuclei compared to the widely used Born approximation. Further progress 
in this direction should help to make experimental results for nuclear effects in DIS more 
reliable. 

Although originally our work has been motivated by the unusual nuclear effects in DIS 
observed in the HERMES experiment |T^| (see Introduction), we did not find anything 
special about the energy and x range of this data. Whatever deviations from the Born 
approximation happen due to the long-range Coulomb forces, they should have similar 
magnitude either at the energy of HERMES, or NMC experiments. 
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Appendix A Calculation of the vector d 



We present here the details of calculations of the vector d in Eq. fl38|). Since we are going to 
use the results also in Sect. |3J] where the amplitude contains two terms in Eq. fl43|) , we will 
treat the photon as a massive particle with mass m 7 which is either zero [the first term in 
(|43|) 1, or equals to my [the first term in fl43|)]. Correspondingly, in all equations of Sect. || 
one should replace v =^ v = ^jv 2 — m 2 . 

The expression for the vector g(pi,p 2 , A) in Eq. (P§|) can be represented as 



+ 



A-niN 



Ptfl V • p 2 ~ Pie 2 V • Pi ) I(q,pi,p 2 ,X) 



q=p 1 -p 2 



Q 2 + + A 2 - 2i& 
( 2iZa(e 1 p 2 - e 2 p{, 



Q 2 + m 2 + A 2 2 
(pi +Pi + v + «A) 2 

$(x) - 2 



• ~ -| iZa 
llV 



ei e 2 , 
D[ - W 2 ) {pm 



{ pi + p 2 + v + i\ 
The following notations are used here, 

$(x) = F(iZa, 1 — iZa; l;x) ; 



p 2 p 1 )(l-£)$'(x) .(A.l) 



(A.2) 
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Sf [X) 

X 


- —<f>(f) ■ 

= X q=p!-P2 = 


^9 A 9 

Q 2 - A 2 


+ + — Liva 


Q 


= q -v ; 


A 2 


= (Pi - P2) 2 




D 1 


— u q=pi-p2 — 


p\-{P2 + v + *A) 2 ; 


D 2 


= Vq=px-p2 — 


( Pl + i> + i\) 2 - p 2 2 . 



(A.3) 
(A.4) 

(A.5) 
(A.6) 

Disregarding the difference between p 12 and ei j2 (equivalent to the dropping off the 



corrections of the order of m jp\ 2 ) we can write, 



2ei 2e 2 



2(z> + iA) 



Di D 2 (pi - p 2 ) 2 + (A - iv) 2 p 1 +p 2 + u + i\ 
Making use of this relation we get, 



(A.7) 



2(A) = 



r Q 2 + m 2 + A 2 -2ii>A 



iZa 



X 



(Pi +P2 + v + iA)2 
'l-5)$'(5) + - Za$(£) + 2 ' (1_£)$,(£) 



Notice that 



and 
where 



in N (eipi - e 2 pi) 
g 2 + m 2 + A 2 - 2iz>A 

2(A-w>) 
Pi -P2) 2 + (A - w) 2 

2(A-«>)(1 -x) 
(pi-p 2 ) 2 + (A + 2z>) 2 ~ g 2 + m 2 + A 2 - 2^ A 

2Za $(5) + 2i(l - 5) $'(5) = i(Za) 2 $(x) , 



2(A — ii/) 



Pi + P2 + ^ + *A 



1 <9x 
x d A 



$(x) = F(l + iZa, 1 - iZa; 2; x) . 
Therefore, the vector g(\) gets the form, 



. (A.8) 

(A.9) 

(A.10) 
(A.ll) 



2(A) = 



AnN(e 1 p 2 - e 2 pi) 

Q 2 



Q 2 + m 2 + A 2 - 2iz>A 
Pi +P2 + v + ^A) 2 



X 



<9$ ^ i(Za) 2 5$(x) 



<9A pi + p 2 + «A + v 



(A.12) 



Here the second term in the brackets of the last factor is small, ~ 9 2 ln(l/#), relative to the 
first term (see below), therefore, it can be neglected for the small scattering angles 
In this case the problem under consideration is reduced to the integral, 



00 

JdXg(X) 



AttN 

— - {e lP2 - e 2 Pi) L 



Q 2 



(A.13) 
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where 



and 



Q 2 + ml + A 2 - 2wA 



Pi + T>2 + i> + «A 
Integrating in Eq. (|A.14[ ) by parts we get, 



L = t iZa $(x 



A=oo 



A=0 



Q j-iZa 



(A.14) 



(A.15) 



(A.16) 



Since x = at A — > oo and $(0) = 1, we have 

A=oo 



A=0 



where 



t = t(A=0) = 

x(A = 0) 



x 



l) iZa - 4 Zq $(xo) 



Q 2 + m 2 

(Pi +P2 + V) 2 

Q 2 

Q 2 + ml 



(A.17) 



(A.18) 



To calculate the rest integrals we notice that the value of x is tiny, except the region of 
A < (Q 2 + m 2 )/(2z>) <C 1. Then, we split the integral in Eq. (|A.16|) into two parts, 



OO Al OO 

J d\ = J dX + J dX , 



(A.19) 



where we chose (Q 2 + m 2 / )/(2z/) <C Ai <C u, so that \x\ <C 1 in the second term in (|A.19|) . 
Therefore, we can fix $(x) = 1 in the second integral in ( |A.19| ) and get, 



a 4.iZa 



+ o(^ 



. 1 yz a _ t iz a + f^ 



(A.20) 



where t\ = £(Ai). 

In the first integral in QA.19 ) we can make use of smallness of A in the expression (|A.15|) 
and rewrite it as, 



t 



Q 2 + m 2 - HXv Xo t 1 



(Pi +P2 + V) 2 

Then the first integral in ( |A.19j ) takes the form, 

Ai 



X X 



(A.21) 



dX-^y$(x) 



iZa (xq to 



§(x) 



X'(| 



x t y Za x~ iZa F(iZa,iZa;l;x) 



t[ Za -ti Za F(iZa,iZa;l;x )+O[ — \ . (A.22) 



Ai' 



Here we used the definition of from Eq. 

Now, taking into account that \x(X = Ai)| <C 1 we eventually arrive at the following 
expression for the integral Eq. ( |A.14| ), 

L = t iZa [F(iZa, iZa; 1; x ) - F(iZa, 1 - iZa; 1; x )} . (A.23) 

Taking into account the relation, 

F(iZa, iZa; 1; x ) = F(iZa, —iZa; 1; x ) + iZa -^(1 + 1 — iZa; 2; x ) , (A. 24) 
we can integrate g(X) as, 



d\g(X) 



Q 2 + ml 



Q 2 + m 2 : 



(Pi +P2 + v) 2 



iZa 



iZa F(l + iZa, 1 - iZa; 2; x ) . (A.25) 



Substituting this expression into Eq. (|38|) we eventually arrive at the final expression for d, 



AttN 



Q 2 + m. 



- \€ X P2 - e 2 pij 



Q 2 + ml 



(Pi +P2 + V 



r,\2 



x [F(iZa, — iZa; 1;xq) — iZa(l — xq) F(l + iZa, 1 — iZa;2;x )] . (A. 26) 



Appendix B Eikonal approximation versus the exact 

calculations 

Here we calculate the total cross section of photoproduction of electron-positron pairs off 



atoms at high energies relying upon the eikonal approximation of ref. [20| and compare 
with the exact results of |2T[ . 



According to [^(J the photoproduction cross section has the factorized form, 



a^Z — > e + e Z) = du d 2 pa(p,u) ^(p,u 



(b.i; 



where a(p, u) is the cross section of interaction of the e + e~ dipole with the atom which 
depends on the transverse separation p and the fraction u = (E + p£)/2uj of the photon 
light-cone momentum carried by the electron. ty(p,u) is the e + e~ wave function of the 
photon, 



am 



{K 2 (mp) + 



u 2 + (l 



u 



K 2 (mp)} 



(B.2) 
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where m is the electron mass, a is the fine structure constant. 
If the atom remains intact, then according to |2"0fl , 

a(p, u) = 2 J d 2 b jl — exp i x(b — up) — i P — up) j 



(B-3) 



where x(b) is defined in 

Replacing in ( B.3|) 6 + (1/2 — w)p =^> 6 we arrive at, 

a(p,u) = 2 d 2 b {l -exp[ix(6 + ) -«%(&_)]} , 



(B.4) 



where 6± = b±p/2. Thus, we conclude that the dipole cross section Eq. ( B.4|) depends only 
on the transverse separation p. 

The phase shifts x(&±) can be expressed in terms of the transverse density of electrons 
in the atom n(s), 

X{b±) = 2Za J dssn(s) In (j-^j , (B.5) 



where b± = \b±\ and the density is normalized as, 

dssn(s) = 1 . 



(B.6) 



The atomic size Rz ~ 1/ {maZ 1 ^) is much larger than the transverse e + e separation 
p ~ 1/m, therefore, one can split the integral in (|B.4[) into two parts, 



OO 2lT 

a{p) =2 J dbb 

n 



1 - exp(«x+ - ix-) = +°2{p) 



(B.7) 



where is the azimuthal angle between b and p; x± — x(b±), and 

b 2vr 



(P) 
0"2 (p) 



2 / 





oo 2n 

2 J dbb 

b o 



1 - exp(ix+ - iX- 
1 - exp(ix+ - %X- 



(B.8) 
(B.9) 



We chose the value of bo satisfying the condition 

— < b < R z . 

m 



(B.10) 



Starting with Eq. ( |B.9| ) for cr 2 (p) we note that the mean value of the e + e separation is 
much smaller than the impact parameter, 



p ~ — ^b <b 
m 
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Therefore, 



X+ ~ X- 



x(b+ -x(b- P ~)=p- V bX (b) +0[ p - 



p COS( 



dx(b) 
db 



2Zacos</> - / dssn(s) 



(B.11) 



Corresp ondingly, 



02 0) 



4vr / 



1- Jo p 



= 4vr(Za)V / 

bo 

Integrating this expression by parts we get 



db 

T 



dx 
db 



ds sn(s) 



+ 0[ 



Ri 



(B.12) 



&0 



db 
T 



oo 




y ds sn(s) 

.b 









ds sn(s) 

- 2 J dbb lnf- J n(&) Jdssn(s). 



(B.13) 



The bottom limits &o of the integrals in the right-hand side (r.h.s.) of ( |B.13| ) can be 
replaced by zero, b =^> 0, with an accuracy of the order of b\jR 2 z <C 1. Then, taking into 
account the normalization condition Eq. ( |B.6| ) we get, 



a 2 (p) = A7t(Za) 2 p 2 



In 



+ (I 



where 



oo / 1 \ °° 

In a = — 2 y db b In ( — | n(6) y ds sn(s) . 



(B.14) 



(B.15) 



Apparently, the value of a ~ i?z depends only on the details of the atomic structure con- 
tained in the electron density distribution n(s). 

Now we turn to the first term 0i(p) in (|B.7| ). Since b < b <C R z , also b± = \b ± p/2\ < 
b <C i?^. Making use of these relations we get, 



X± 



2Za / dssn(s) In 



b± 
oo 

2Za y dssn(s) In 



it!! 



(B.16) 
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Consequently, 



X+ ~ X- 



2Za In — dssn(s) + Ol-£ 



Za In 



b + p/2) 2 



Za In 



b 2 — bpcoscj) + p 2 /4\ 
b 2 + bp cos (f) + p 2 /4 J 



(B.17) 



Substituting this expression in (|B.9|) we get <Ji(p) in the following form, 



60 2tt 

cxi(p) = 2 /g?66 



1 



b 2 — bpcoscj) + p 2 /AY Za 
b 2 + bp cos <fr + p 2 /4J 



(B.18) 





To proceed to the further modifications of this expression we employ the relation, 

1 -X Y = -XYF(1-Y, 1;2;1-X) , (B.19) 

with 

b 2 — bp cos <f) + p 2 /4 



b 2 + bp cos (ft + p 2 /4 
Y = iZa . 

We also make use of the following representation for the hypergeometric function, 



F(l — Y, 1; 2; 1 — X) 



I dt ' {l - tY 

r(i-y)r(i + y)7 i-t(i-x) 

Applying these relations to (B.18| ) we get, 

1 

2 



(B.20) 



(B.21) 



T(l - iZa) T(l + iZa) 



dtr lZa (i-ty Za L(b ,p,t) 



iZa 



(B.22) 



where 



L(b ,p,t) 



b 2vr 
dbb 

n b 



2b p cos < 



b 2 + (1 - 2t)6pcos</> + p 2 /4 



2vrp 2 (l - 2t) 



x 



. b 2 + p 2 /4 + VD 2 



In 



ft 2 + [1 - 2(1 - 2t) 2 ]p 2 /4 + VP 
2p 2 t{\ -t) 



and 



D 



b 2 + P 



2\ 2 



4 



6 2 p 2 (1 - 2ty 



(B.23) 



(B.24) 



According to the above convention p <^ bo and this expression can be further simplified, 



L(bo,p,t) = 7ip 2 (l-2t) 1-ln 



h 2 



P 2 t(l-t) 



O 



Pi 

bl. 



(B.25) 



22 



Now integration in ( B.22 ) can be performed analytically, and we arrive at the final 
expression for <7i, 



(7i (p) = An(Za) 2 p 2 



In 



V 



- 2 -Re^(l + iZa) - C 



where 



= — lnlYiu) 
aw 



(B.27) 



C 



-V(l) = 0.5772 . 



Adding Eqs. ( p.!4j) and (p. 27| ) we eventually get the dipole cross section ( P-12|) in the 
form, 

a(p) = o Born (p) - Acr(p) . 



Here 



0" Born (p) 



4tt(Z«)V 



In a - 



and 



4vr(Za)V Re^(l + iZa) + C = 47r(Za) 2 p 2 f(Za 



f[Za) = {Zaf £ wrim 



(B.28) 

(B.29) 
(B.30) 

(B.31) 



Now we are in position to calculate the photoproduction cross section substituting 
Eq. qB^D into ®, 



cr(7Z — > e + e Z) = <7Born(lZ — > e + e Z) — Act(7Z — > e + e Z) 



where 



a jBor . n (7Z e + e Z) 
Aa(7Z -> e + e"Z) 



28 Z 2 a 3 
9 m 2 

28 Z 2 a 3 
9 m 2 



ln(2am) - C 
f(Za) . 



83 
42 



(B.32) 

(B.33) 
(B.34) 



Note that the Coulomb correction Aa^Z — > e + e~Z) in Eq. ( |B.34| ) coincides with one 
calculated in [fjlj using a very different technique. At the same time, the Born terms are 
different what one should have expected. Indeed, the Bethe-Maximon theory describes pho- 
toproduction of e + e~ pairs off a point like nonscreened nucleus, while we are dealing with pair 
production off atoms. Equivalence of the Coulomb corrections calculated within the eikonal 
approximation and in the Bethe-Maximon theory in the case of photon bremsstrahlung can 
be proven exactly in the same way, since both processes are controlled by the same dipole 
cross section. 
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